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We investigate the relation between the E-, spectrum in B — > X s ^y decay and the P+ spectrum in 
semileptonic B —* X u £i> decay (P+ is the hadronic energy minus the absolute value of the hadronic 
three-momentum), which provides in principle the theoretically simplest determination of \Vub\ from 
any of the "shape function regions" of B — » X u lv spectra. We calculate analytically the P+ spectrum 
to order ati(3o, and study its relation to the B —* X a ^ photon spectrum to eliminate the leading 
dependence on nonperturbative effects. We compare the result of fixed order perturbation theory 
to the next-to-leading log renormalization group improved calculation, and argue that fixed order 
perturbation theory is likely to be a more appropriate expansion. Implications for the perturbative 
uncertainties in the determination of \V u b\ from the P+ spectrum are discussed. 



I. INTRODUCTION 

The determination of the magnitude of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element V u b via in- 
clusive decays is theoretically involved, because the ex- 
perimental cuts required to suppress the B — > X c lv back- 
ground tend to restrict the B — > X u iD phase space in a 
way that gives rise to theoretical complications [1]. The 
local operator product expansion (OPE) [2, 3] breaks 
down in the regions Eg > (m% — m 2 D ) / (2niB) and 
nix < tud [4, 5] (where Eg is the charged lepton energy 
and rrix the hadronic invariant mass), because numeri- 
cally m 2 ~ Aqcd?ti&- In these regions an expansion in 
in terms of nonlocal matrix elements is still 
possible [6], and the leading term can be measured in 
B -> A s 7 decay [7]. (The local OPE is valid in the 
q 2 > (m B — m D ) 2 region [8], q 2 being the lepton invari- 
ant mass, but there are other issues for that cut [1].) The 
kinematic region, in which the final hadronic state has 
high energy ~ nib but low invariant mass ~ \/AqcdW6) 
is typically known as the "shape function region." 

It was pointed out [9] (see also [10, 11]) that sep- 
arating b — > u from b — > c using another variable, 
P+ = Ex — \px\, where Ex and px are the energy 
and three-momenta of the hadronic final state, may pro- 
vide advantages compared to Eg or nix- At lowest 
order in AQCD/^b dr B ^x u iu/dP+ 1S proportional to 
drs_ > x S 7/d-B 7 evaluated at E 1 — (tub — -P+)/2 (since 
P + in B — > X s j equals mg — 2E 1 ). Thus, to predict the 
B — > X U £D rate in the P + < m 2 D /niB region that is free 
from the charm background, we only need to know the 
B — > X s 7 photon spectrum for E 1 > (m 2 B — m 2 D )/ (2dib); 
i.e., a 330 MeV region near the endpoint, which is already 
precisely measured [12-14]. This is also what is needed 
to predict the Eg endpoint spectrum, but it is signifi- 



cantly smaller than what is required to determine the 
mx spectrum [15, 16]. 

A convenient way to express the relation between the 
P+ spectrum in semileptonic b — > u decay and the photon 
energy spectrum in B — * X s ^f decay is to relate weighted 
integrals of the two spectra [17] (see also [10]). One can 
write 

where we have defined P 7 = nis — 2_B 7 , and at lead- 
ing order in Aqce) / m b the weight function W(A,P~) is 
calculable in perturbation theory. 

In the shape function region, A ~ Aqcd, the pertur- 
bative expansion of W contains logarithms of the ratio 
of scales y / mbA~QCD /???,&. For nib ^> Aqcd, these loga- 
rithms are large and can spoil the convergence of per- 
turbation theory, so must be resummed using renormal- 
ization group techniques. This can be done using tra- 
ditional perturbative QCD methods [17-20] or by using 
the Soft-Collinear Effective Theory (SCET) [21]. Cur- 
rently W can be extracted from known results to leading 
log (LL) and next-to-leading log (NLL) accuracy. How- 
ever, for the true value y/ nib Aqcd / mb ~ 1/3, it is not 
clear that the leading log expansion is appropriate, and 
a fixed order calculation in a s may give a better approx- 
imation to W } 

In this paper we address this issue by calculating the 
order ev 2 /3o corrections to the P+ spectrum, where /3q = 



This is reminiscent of the resummation of logs of m c /mi, ~ 1/3 
in exclusive B — > D*lu decay at zero recoil, in which the leading 
log calculation is a poor approximation to the one- or two-loop 
results [22, 23]. 
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11 — 2n//3 is the first coefficient of the QCD /3-function 
(we shall refer to this order as BLM [24] ) , from which we 
determine the corresponding corrections to the weighting 
function W. Since f3o ~ 9 is a large number, such terms 
dominate the two-loop corrections to many processes [24, 
25]. We find that a considerable part of the two-loop 
expression arises from terms which in the renormalization 
group improved perturbation theory only occur at next- 
to-next-to leading log (NNLL) order, and we argue that 
fixed order perturbation theory is a more appropriate 
expansion. We discuss the implications of this for the 
uncertainty in the determination of \V ub \. 

II. THE SPECTRA AT ORDER a^f3 

We first present analytical results for the parton level 
spectra in B — > X u lv and B — ► X s j to order a 2 f3o- 
Our results for the P + spectrum are new, whereas those 
for the photon spectrum arc collected from the existing 



literature for completeness. 

At the parton level, the appropriate variables are 

p+ = (v-q/m b )-n, x=l-2E 1 /m b , (2) 

for B — > X u lv and for B — ► X s -y decay respectively, 
where n is a light-like four- vector in the direction of —q, 
v is the four- velocity of the B meson. These variables are 
simply related to the experimentally observed hadronic 
variables by 

P+ = E x - \px | = (m B v -q)-n = m b p+ + A , 

P 7 = m B — 2_E 7 = m b x + A , (3) 

where A = tub — m b . 



A. The p+ spectrum in B — > X u lu 

The p+ spectrum to order a s is given by [9] 
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where T = G^|V^t| 2 m b /(1927T 3 ), Cp — 4/3, and m b is the b quark pole mass. The * distributions (for n > integers) 
arc defined by 
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Integrating Eq. (4) over < p+ < 1 reproduces the total rate 1 — (o^Cf / (47r))(27r 2 — 25/2). 

The BLM part of the two-loop result may be obtained from the one-loop calculation with an arbitrary gluon 
mass using the method of Ref. [26]. We calculated the spectrum for p+ ^ 0, for which only bremsstrahlung graphs 
contribute, and then determined the coefficient of S(p + ) by comparing with the total rate [27]. The result is 
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where Li 2 (z) = — j Q z dfln(l — t)/t is the dilogarithm, and 
Li 3 (z) = J*dtU 2 (t)/t. 

B. The x spectrum in B — > X s 7 

We concentrate on the part of the spectrum that arises 
from the operator Oj = (e/167r 2 ) mt sl u^F^ &r. This 
gives rise to the dominant part of the photon spectrum, 
and while other operators also influence the spectrum, 



the photon spectrum is only known analytically to order 
a 2 (3o for this piece (and for 8 , but not for the terms 
involving the four-quark operators). To study the con- 
vergence of the perturbative expansions and assess the 
theoretical uncertainties in the P + spectrum, it is suffi- 
cient to consider O7. However, ultimately, for the actual 
determination of \V u b\ the other contributions should also 
be included. 

The photon spectrum to order a s is given by [28, 29] 
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where T 7 = G 2 F | V tb V t * s | 2 a om mf [m b (m b )CY S (m b )] 2 / (32tt 4 ) , and we have set (i = m b for convenience. The BLM 
correction to the photon spectrum due to O7 may be obtained by combining results given in [30] and [31]; again, 
setting fi = m b it reads 
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III. RELATION BETWEEN THE SPECTRA 

One can define a weighting function, W, that relates, 
in the shape function regions, weighted integrals of the 
photon energy spectrum in B — > X s j to integrals of the 
P+ spectrum in B — ► X u lv: 
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The weighting function W can be computed perturba- 
tively at leading order in Aqcd/w{,, because the shape 
function contribution that captures the leading nonper- 
turbative physics at scales of order Aqcd is identical in 
the two spectra. The corresponding perturbation series 
depends on the scales m b and ^ym b A^c^>. The function 
W is also free of logarithms of the form a™ ln m (m 2 /yu 2 ), 
m = n + 1, . . . , 2n, which are universal in the two spectra 
and cancel from the relation in Eq. (9). 

The origin of the m^/m 2 factor in the definition of 
W deserves comment. At leading order in Aqcd/"J6, the 
effects of the shape function f(u>) may be simply included 
by making the replacement 



m b 



ml = m b + u , 



(10) 



differential rate with the shape function f(u>) [7], 

ar = / <ir !M " ' ;; | ^ . /m<U'. (11) 



= [ dr parton | . f(w)dw. 



This prescription also generalizes to higher orders in 
Aqcd/w;,, provided that the replacement (10) is only ap- 
plied to factors of m b which arise from kinematics, and 
not from coefficients of operators in the Hamiltonian [16] 
and correctly reproduces the class of "kinematic" sub- 
leading effects (proportional to the leading order shape 
function) which do not arise from the expansion of the 
heavy quark fields. 

Applying this procedure to the tree level spectra, we 
have 



dT„ 
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in the tree level partonic rate, and then convoluting the 
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cx (m b + to) 5 5{P+ -A + uj) f(oj) duj 
'+ J 

= J(m B - w) 5 5(P+ - u>) /(A - w) d£ (12) 
for semileptonic decays (where uj = A — u>), and 

dr f 

cx m b (m b ) 2 / (m B - w) 3 5(P 1 - u>) /(A - uj) dw 

(13) 

for radiative decays, where two powers of m b (m b ) origi- 
nate from C 7 in the effective Hamiltonian to which the 
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replacement (10) does not apply. This gives 
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Thus, it is natural to pull the leading factor of 
m 2 B /m b (m b ) 2 out of the definition of W. As we will show, 
factoring out this term rather than leaving the "partonic" 
factor m 2 /raj(mt) 2 in the definition of W dramatically 
reduces the absolute size of the perturbative corrections 
to W. 



To calculate W, we first expand both parton level spec- 
tra in powers of p + and x respectively, since in the shape 
function region both are 0(Aqcd/w£,). At leading or- 
der in Aqcd /m b , this simply corresponds to keeping the 
first four terms in Eqs. (6) and (8). At leading order in 
the Aqcd jmb expansion relevant for the shape function 
regime, the Feynman diagrams which give the coefficient 
of the shape function give the partonic result with the 
substitution p + — > p + + k + and x — > x + where k + 
is the light-cone component of the residual momentum 
of the heavy quark, which is of the same order. (The 
dimcnsionless variables k+ and <2> have support between 
— oo and A/m b .) Thus, we have 
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and similarly for dl^/dx. This gives, for each spectrum, the matching conditions onto the leading nonlocal operator 
Oq{uj) = bS(uj — iD + ) b in the OPE. The B meson matrix element of Oo(u) gives the leading order b quark light-cone 
distribution function, or shape function. In principle, to determine the coefficient of Oo(^) we must include the 
radiative corrections to the parton level matrix element of Oo(u>); however, since these terms are common to both 
spectra and therefore drop out of W, we do not need to worry about them here. 

Including the full one-loop corrections and the BLM two-loop contributions, and setting the scale \i = m b for 
simplicity, the function W is given by 
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where the ellipses denote non-BLM two-loop terms, higher orders in perturbation theory and nonperturbative correc- 
tions suppressed by powers of AQCD/m b . Eq. (16) is the main result of this paper. 

It is instructive to compare Eq. (16) with the next-to-leading log result in SCET. Using the results of Rcfs. [10, 32] 
(see also [33]), we find 
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ip' is the derivative of the digamma function, tjj(z) = T'{z)/T{z), and 
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Here Tq = AC F and Ti = 8CV(67/6 — 7r 2 /2 — 5n//9) are the first two coefficients of the cusp anomalous dimension, 
and Pi — 102 — 38n//3. Setting ^ ~ AQCD^-b sums all leading and subleading logarithms of the form a™ hi 1 {iii / m b ) 
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and a™ In™ (fj,i/m b ). Expanding Eq. (17) to order a 2 (m b ), we find 
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where the first line is the leading log result, and the second and third lines are subleading logs. Note that the 
renormalization group equations sum logs of m 2 /^i 2 , not m b (A — P 7 )/fj, 2 , since for /z 2 ~ AqcdTO6 only the first are 
parametrically enhanced. The result is formally independent of fa, as can be seen by re-expanding the logs in Eq. (20), 
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where we have dropped terms of order a 2 \n n [m b {A — P 7 )/^f], which are next-to-next-to- leading order. This result 
agrees with the corresponding one-loop and two-loop BLM terms in Eq. (16). 

Equations (16) and (21) both provide approximations to the full expression for W at two-loops, so it is instructive 
to compare them. Numerically, the 0{a 2 ) terms are 
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where a complete two-loop calculation is required to de- 
termine Co- 

For both the double and single log terms, the BLM- 
enhanced term is about a factor of two larger than the 
non-BLM term, which suggests that it may also dominate 
the nonlogarithmic term to a similar degree. In contrast, 
the leading log approximation is clearly poorly behaved 
at this order: for ^ 2 /m 2 <~ (A — Py)/m b = 1/9, the 
double, single and non-logarithmic terms in Eq. (22) are 
in the ratio 

©(log 2 ) : O(log) : O(log ) = 1 : 0.87 : (-0.86 - 0.02c ) . 

(23) 

This reflects the fact that the logarithmic enhancement 
is not sufficient for the double log to dominate over the 
single or zero log terms, nor for the single log to dominate 
over the BLM-enhanced piece of the zero log term. The 
same conclusion is reached by comparing the size of the 
LL and NLL terms in Eq. (20). This suggests that a fixed 
order calculation, rather than a leading log calculation, 
is more appropriate for W. We will discuss the numerical 
significance of this for the extraction of \V ub \, along with 
the poor convergence of both the leading log and fixed 
order calculations, in the next section. 

Another source of corrections to W comes from terms 
suppressed by AQCD/'Tifc- Some of these modify W with- 
out introducing new unknown hadronic matrix elements, 



while others involve nonperturbative matrix elements of 
higher-dimension nonlocal operators [6], which cannot 
presently be computed. The mismatch of the powers of 
m b in the two spectra gives rise to the order AQcr>/m b 
correction to W contained in Eq. (14), — 2P 7 /m_B- At the 
same order, there are also corrections that come from ex- 
panding the b quark fields in powers of Aqcd /m b . These 
effects are sensitive the Dirac structure of the current, so 
do not cancel from W. We can extract the full Aqcd/^6 
correction from Ref. [34], and find 

w(AqcdH = 8P 7 - 2A 

3m B '"' v ' 

The corrections that depend on subleading shape func- 
tions are also given in [34] and involve five unknown non- 
perturbative functions, whose effects we do not attempt 
to model here. Some of the 0(a s ) corrections to Eq. (24) 
can be obtained by expanding the results in Sec. II to 
higher order in p + and x. We find that these are small 
compared to the perturbative uncertainties in W at lead- 
ing order in AQcv/m b . 

IV. IMPLICATIONS AND DISCUSSIONS 

To eliminate the charm background kincmatically from 
semileptonic B — > X u iv decays, one has to impose a cut 
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FIG. 1: W(A,P 7 ) as a function of P 7 . At tree level W = 1, and the curves include different orders in the expansion: (a) the 
0(a s ) result; (b) all known 0(a 2 ) terms, Eq. (22) with co = 0; (c) the LL RGE resummed result; (d) the NLL RGE resummed 
result, Eq. (17); (e) the C(a^/3 ) result, Eq. (16) and (f) the NLL expression expanded to 0{a 2 3 ), Eq. (20). 



P+ < mj-f/rriB — 0.66 GeV. Then measurements of the 
P+ spectrum in B — > X u lv and the P 7 spectrum in B — > 
X s j, together with the theoretical input of W(A, P 7 ) will 
allow a determination of \V u b\ using Eq. (9). 

In Fig. 1 we plot IT(A,P 7 ) for A = 0.66 GeV, 
a s (rrib) = 0.22 and nib = 4.8 GeV in different approxi- 
mations. At tree level, W(A, P 7 ) = 1. Curve (a) is the 
order a s result, while the result with all known a 2 , contri- 
butions is shown in curve (b) [i.e., Eq. (22), the BLM re- 
sult plus the non-BLM coefficient of the double and single 
logs] . The dashed black curves show the RGE resummed 
results for fa = \Jrrib A ~ 1.78 GeV. The straight line (c) 
is the LL result, whereas the long dashed curve (d) is the 
NLL result in Eq. (17). The LL and NLL results use 1- 
and 2-loop running for a s , respectively. Finally, the dot- 
dashed curve (e) contains just the BLM terms at order 
a.g, Eq. (16), while the short dashed curve (f) contains 
NLL expression expanded to order a?, Eq. (20). 

The difference between the dashed (d) and (f) curves 
shows that the terms in the NLL sum beyond 0(a 2 s ) are 
not negligible. However, the logs that are summed do 
not dominate over other contributions in the perturba- 
tion series; e.g., the difference between Eqs. (20) and (21) 
is NNLL, but it is comparable to the non-BLM LL and 
NLL terms at G(a 2 s ). Therefore, we view the fixed order 
result, curve (b), as our best estimate of the perturbation 
theory prediction of W. The difference of these curves 
provides an estimate of the uncertainty related to higher 
order uncalculated terms. In addition, the /Lti-dependence 
of the NLL resummed result is larger using two-loop than 
one-loop running for a s , also indicating that renormal- 
ization group improved perturbation theory may not lead 
to an improved expansion. 

To assess the significance of these results for the 
determination of \V u b\ using Eq. (9), we integrated 
these various expansions of W(A, P 7 ) against a simple 



TABLE I: Weighted integral J (| A dP 7 W(A, P 7 ) (dr s /dP 7 ), 

normalized to f Q A dP 7 (dI\,/dP 7 ) for A = 0.66 GeV, taking 
the simple parametrization (25) of the experimental photon 
energy spectrum. 
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parametrization of the experimentally measured B — > 
AT S 7 spectrum [13]. For the optimal cut, A = 0.66 GeV, 
the integral on the right-hand side of Eq. (9) normalized 
to that integral at tree level (corresponding to W = 1) 
is shown in Table I at order a s , ce^fio, using the LL and 
NLL RGE resummations, and all known 0(a 2 ) terms. 
The simple parametrization 



dP, 
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_2 
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with A ~ 0.9 GeV and a ~ 6.1 provides a crude but, 
for our purposes, sufficient fit to the data, 2 as the ratios 
in Table I are quite insensitive to the precise shape of 
the spectrum. (Using the parameters A ~ 0.66 GeV and 
a ~ 3.3, which is the shape function fit rather than the 
photon spectrum [35], and so corresponds to a rather 



2 A parameterization of the physical B — > X s -/ spectrum is not 
available, only the raw spectrum or fits to shape function pa- 
rameters [35]. However, the shape function fit includes 0(a s ) 
corrections differently than how they enter W. Thus, the num- 
bers in Table I should be taken only as indicative of the size of 
the corrections. 
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different shape, only changes the entries in Table I to 
1, 1.08, 1.15, 1.10, 1.18, 1.14, respectively. However, it 
changes J Q A dP 7 (dT s /dP^) by about a factor of two.) 

The convergence of the result in Table I going from tree 
level to 0(a s ) to 0(a 2 /3o) is poor, and that going from 
tree level to LL to NLL resummation is worse. This poor 
behavior of the perturbation series may be related to the 
fact that in the nonlocal OPE there are 0(Aqcd/™&) 
nonpcrturbative corrections to Eq. (9) arising both from 
the explicit factor of A in Eq. (24) as well as the sub- 
leading shape functions, and so there is a renormalon 
ambiguity in the perturbation series for W at this order. 

The subleading twist terms in Eq. (24) give a —0.19 
correction to the numbers in Table 1 . (The "trivial" part 
of this correction, — 2P 7 /ms, that is due to the mismatch 
of the two powers of m& smeared in the two spectra, gives 
numerically the same result.) While this is a calcula- 
ble effect and not an uncertainty, there are incalculable 
subleading shape functions that enter at the same or- 
der, which can only be modelled. Thus, the perturbative 
and nonperturbative uncertainties in this determination 
of | V u b | are probably comparable. 

Because the experimental photon spectrum is peaked 
around P 7 ~ 0.8 GeV, the weighted integral (9) is dom- 
inated by small values of A — P 7 , increasing the im- 
portance of the logarithmically enhanced terms in (21). 
However, this numerically large logarithm is not summed 
by the RGE, since it is not logarithms of (A — P 7 )/m6 
but rather of fii/nib which are summed. 

Had we not factored out m^/m 2 in the definition of 
W in Eq. (9), there would be an additional perturbative 
factor of m^/m 2 in W which, when expanded out, would 
modify the expressions in Eqs. (16), H(a) in (18), (20), 
(21), (22), and (24). If we define a weighting function 
W'(A, P 7 ) in this way, then instead of Eq. (22) we obtain 



C F a 2 s {m b ) 
(4tt) 2 



(0.83A, + 3.41) In 2 - m * (26) 



A-Pv 



+ (4.67/3 - 1-35) In 



TOfc 



A-P, 



+ (32.3/30 + ^) 



and the non-logarithmic BLM term is considerably larger 
than cither of the logarithmic terms. In this case the 
numerical results in Table I would read 1, 1.29, 1.51, 
1.10, 1.43, 1.53, respectively, and we would have to 
assign a much larger perturbative uncertainty to W 
than to W . However, this is due to the bad perturba- 
tive behavior of m 2 /m 2 , and not of the spectra them- 
selves. In this case the analog of Eq. (24), V^ /(a « cd/ " i) = 
— 8(P 7 — A)/(3ms) + . . ., gives a —0.01 correction to the 
figures in the previous sentence. It is interesting to check 
the consistency of the results. Combining all known order 
a 2 and Aqcd/w& terms, we obtain 



2 f A dP 7 W(A, P 7 ) (dr 77 /dP 7 ) 

J A dP 7 H/'(A,p 7 )(dr 77 /dP 7 ) 



(4.27 GeV) 2 , 



(27) 



quite consistently with the physical value of m^mj). In 
comparison, Eq. (27) with the NLL result for W and W' 
gives (4.50 GeV) 2 . We learn that if we keep all two-loop 
corrections, the physical result is quite independent of 
whether we calculate it in terms of W or W, while the 
same cannot be said about the NLL resummation result. 

The measured B — > X s j photon spectrum together 
with W(A,P 7 ) given by the sum of Eqs. (22) and (24) 

determines / Q A dP 7 W(A, P 7 ) (dP s /dP 7 ), which in turn 
determines \V u b\ from a measurement of the partially in- 
tegrated P + spectrum in B — ► X u iv using Eq. (9). The 
theoretical uncertainty of \V u b\ from perturbation theory 
alone using this method is half the error of the results 
in Table I. However, because of the poor behavior of 
the perturbation series, the full two-loop calculation of 
VF(A, P 7 ) would be desirable. Furthermore, the pertur- 
bative series is likely to improve if the unknown matrix 
elements at C(Aqcd/to;,) are expressed in terms of phys- 
ical quantities, so that the leading renormalon ambiguity 
in W is cancelled. In addition, and probably more im- 
portantly, effects of operators other than 7 need to be 
included, in particular that of O2 and Os may be im- 
portant. In the NLL resummed result, including these is 
straightforward using Eq. (4) in [32]. However, for the 
BLM result one needs to combine the analytically known 
virtual contributions [30] with the bremsstrahlung con- 
tributions [31], which are only known numerically for the 
O2 operator. Work in this direction is in progress. 

In summary, we calculated the order a 2 /3o corrections 
to the P+ spectrum in B — > X u iv decay and studied the 
uncertainties in extracting |V^b| using a measurement of 
the P 7 spectrum in B — > X s j. We showed that the 
factor of ml /ml in W(A, P 7 ) at lowest order naturally 
becomes to^/to 2 when subleading effects are included, 
and results in much reduced perturbative corrections. We 
found that the NLL RGE resummation is of limited use, 
because the logs that it sums do not dominate over the 
non-logarithmic terms. This may have implications for 
the phcnomcnological usefulness of other applications of 
RGE resummations in inclusive heavy to light decays. 
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